Mott-Insulator to Liquid Transition and Population Trapping in Ultracold Fermi 
Gases by Non-Equilibrium Modulation of the Optical Lattice 
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An ultracold gas of interacting fermionic atoms in a three dimensional optical lattice is considered, 
where the lattice potential strength is periodically modulated. This non-equilibrium system is non- 
perturbatively described by means of a Keldysh-Floquet-Green's function approach employing a 
generalized dynamical mean field theory (DMFT). Strong repulsive interactions between different 
atoms lead to a Mott-Insulator state for the equilibrium system, but the additional external driving 
yields a non-equilibrium density of Floquet-states and a transition to a liquid or conducting state. 
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In condensed matter physics quantum dynamics can 
only be analyzed with the premise of severe influences due 
to strong coupling to the environment in the experimen- 
tal situation. This aspect makes it challenging to prepare 
bare quantum states far-from-equilibrium in a controlled 
way. Thus ultracold gases are wonderful systems to study 
the pure influences of quantum non-equilibrium effects in 
first instance, since the properties of such systems can 
be tuned almost without restrictions. At this point the 
reader might ask the question "Why non-equilibrium ?" . 
Whereas calculations in thermodynamically equilibrated 
systems are well established in theory and rather com- 
phcated structures requiring high-end numerics can be 
solved, it is still a challenge to determine results for 
driven systems [ll]. Nevertheless those systems are ex- 
perimentally and technologically interesting [2, 3] . Even 
if a system is in the steady-state regime, it still not re- 
sides in thermodynamical equilibrium if it is driven e^. 
by a nonzero current or a difference in the potential y, 
and the reflection of that fact in theory is highly desir- 
able. Besides that aspect it is a prestigious and intriguing 
aim to investigate fascinating quantum properties such as 
quantum entanglement and the onset of decoherence and 
dephasing [5-12| which have been measured experimen- 
tally with spectacular results. 

Recently, intriguing experiments on non-equilibrium 
dynamics of ultracold gases, both fermionic as well as 
bosonic, [l3l-[l8| have become possible. Among many as- 
pects, also the long time limit of such non-equilibrium 
systems has been studied 19|-21|, and applications using 
the non-equilibrium processes as a tool for transport in 
so called quantum ratchets have been investigated [22 1. 

The Hubbard Hamiltonian is one of the most relevant 
models for investigating strongly correlated systems in 
condensed matter theory, of both, bosonic as well as 
fermionic nature (23l-[25|. In this Letter we study the 
characteristics of an interacting ultracold Fermi gas ex- 
posed to a periodic potential modulation of the optical 
lattice with a constant amplitude. This configuration 
corresponds to a stationary non-equilibrium state, which 



requires suitable techniques such as the Keldysh formal- 
ism. The periodic modulation leads to a 'dressing' of 
atoms well known from the application of light fields in 
quantum optics. We discuss a dynamical mean field the- 
ory (DMFT) solution for the Floquet-Keldysh [26] ap- 
proach. The non-equilibrium caused dressed states arise 
as Floquet side-bands. In the Mott-Hubbard-gap we de- 
rive a complicated modulation induced structure of many 
particle states, and therefore a transition from the Mott 
insulting regime to a liquid phase which leads to a fi- 
nite conductivity. The occupation number for these gap 
states is investigated and we find a trapping of population 
which results in an 'inversion' for increasing modulation 
frequencies. 

The Fermi gas in the modulated optical potential is 
schematically shown in Fig. [H The Mott-insulator state 
in equilibrium is characterized by a tight binding model 
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FIG. 1: (Color online) (a) The strong local interaction U 
between atoms in the equilibrium system, causes the single 
energy band to split into the lower and the upper Hubbard 
band, separated by an excitation gap. (b) Ultracold fermionic 
atoms in a 3 dimensional optical lattice (black). The opti- 
cal lattice is periodically modulated with the frequency Ql, 
yielding a minimum (red) and a maximum (blue) potential 
modulation. The periodic modulation introduces an addi- 
tional, time dependent contribution to the local atom energy 
E{t) and the atomic tunneling amplitude T(r), see text. The 
sketched potentials of this non-equilibrium system, therefore 
represent three different snapshots in time of the external pe- 
riodic modulation. 
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with strong repulsive onsite interactions U which result 
in a band splitting and the establishment of the char- 
acteristic Mott-Hubbard-gap. The periodic modulation 
of the optical lattice potential influences the behavior of 
tunneling t from lattice site to the nearest neighbor as 
well as the onsite energies eo at each lattice site, because 
these quantities depend on the lattice potential V{r). 
The equilibrium tunneling amplitude t must be replaced 
by t + T(r), as well as the equilibrium onsite energy has 
to be replaced by eo + E{r). These time r dependent 
terms are periodic themselves. Therefore, we consider 
the following Fermi-Hubbard-Hamiltonian 

Hir) = ^[eo + i?(r)]ct,c,, (1) 

The Index i labels the lattice site and a the spin, {ij) 
implies summation over nearest neighbores, c] ^ and c- ^ 
create (annihilate) a fermionic atom with spin a at lat- 
tice site i. The time dependent contributions due to the 
periodic modulation of the potential are given by 

E{r) = E cosiQLT) T{r) = T cos{nLr) (2) 

where I^l is the frequency of the lattice modulation, r is 
the system time and E and T are the respective ampli- 
tudes of energy and hopping or tunneling contribution. 

The appropriate technique to describe this periodically 
driven system far off thermodynamical equilibrium, com- 
bines both the Keldysh ^ and the Floquet ^ Green's 
function approach. Due to the non-equilibrium charac- 
ter the Green's function truly depends on two time argu- 
ments, and the following two-time Fourier transform is 
employed 

G-^Jk,c^) = /dr.ez:= / dr,^e^(-^^-)-- 

J-oo I J-T/2 
Xe^("'-")"-^-Gf (k,T,ei,T,„), (3) 
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FIG. 2: (Color online) The local density of states (LDOS) 
given by Im ^^^(cj) for the non-equilibrium systems char- 
acterized by T = 2.0,3.0,4.0 (from top to bottom), for 
U = 4.0 and E ^ 1.0. The LDOS is displayed as a func- 
tion of atomic energies uj and lattice modulation frequencies 
Ql, D is the half bandwidth. 



where m, n are the Floquet-indices labeling the Floquet- 
modes of the system, which are interpreted as quantized 
lattice oscillations, the phonons. Therefore the system is 
constrained to absorb and emit energy only in multiples 
of energy quanta H^l- The Keldysh indices are a, /3 = ± 
indicating the branch of the Keldysh contour and T = 
is the time period. The system time is shifted to a center 
of motion time Tcm + ^^^^^ and a relative time coordinate 
^rei = ^1 — For Completeness and later comparison, 
we note that for the non- interacting case, i.e. U = the 
Hamiltonian Eq. ([T]) can be solved analytically, yielding 
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We note that Ck is the dispersion induced by the stan- 
dard hopping t, and ek originates from the modulation 
induced hopping contribution T(r), p is the integer sum- 
mation index. To solve the full, i.e. driven and inter- 
acting system {U 7^ 0) at zero temperature and half 
filling, we generalize a non-equilibrium dynamical mean 
field theory (DMFT) [Uj. The DMFT maps the in- 
teracting system onto a local impurity model [29|, Isoj . 
which is solved here by iterated perturbation theory 
(IPT) technique [3l|, which is a diagrammatic method 
and has also been generalized for non-equilibrium sys- 
tems. This substantial numerical effort results in the 
full knowledge of the non-equilibrium Floquet-Keldysh- 
Green's function, revealing e.g. the local density of states 
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FIG. 3: (Color online) Shown is the cj and k integrated imag- 
inary part of retarded component of the calculated Floquet- 
Keldysh Green's function, i.e. / dcj / d^k^^lm G'o,n(c(;,k) 
as a function of for different strengths of the induced 
hopping T as shown, for energy E — 1 . The onsite re- 
pulsion strength is U = 4.0, cf. Fig. [2] Due to numer- 
ical restrictions, the number of Floquet modes included in 
the actual calculations is limited. Throughout this letter we 
use n = — 10 . . . + 10, i.e. we include the lowest 21 Floquet 
modes in our calculations. This finite number introduces an 
error to the numerical results as ^ 0. / dcjImC^n must 
be equal tt if the solution is valid. The deviations show the 
range of validity. For small lattice modulation frequencies Q l , 
the numerical quality depends on the strength of the induced 
hopping T for a given E. For strong correlations of about 
U — 4.0, i.e. the corresponding equilibrium system is in the 
insulating regime, the deviation is less then 10% for frequen- 
cies Q.L > 0.25 which marks the prospected range of validity 
for the parameters used in this Letter. 

(LDOS), the non-equilibrium distribution function, and 
the relaxation times by means of the self-energy. A solu- 
tion for the LDOS is shown in Fig. [21 where we display 
ImG^(a;) = / d^kX;^^ImG^^(cj,k). We need to ac- 
count for all the emission and absorption events of exter- 
nally induced energy quanta HVLl^ and therefore we have 
to sum over all indices m, n. We set = 1. Furthermore 
all energies /7, eo, T, uo and Q.l are given in units 
of of half bandwidth D. 

In Fig. [2] we discuss the development of pronounced 
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FIG. 4: (Color online) Displayed is the non-equilibrium dis- 
tribution function of the fermionic ultracold gas as a function 
of atomic energy uj and external modulation frequency Ql- 
The parameters for this original insulating state are [/ = 4.0 
,T = 2.0£^=1.0. Please note the behavior around external 
modulation frequencies of — 1, also see text. 
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FIG. 5: (Color online) Displayed is the non-equilibrium distri- 
bution function of the fermionic gas together with the corre- 
sponding imaginary part of the retarded component of the lo- 
cal Green's function (DOS) for parameters U — 4.0 , T — 2.0 
^ = 1.0 at a series of different external modulation frequen- 
cies , G {0.3,0.5,0.7,1.0,1.3,1.5}, from top to bottom. 
The band gap in equilibrium has the width 2D. For a de- 
tailed discussion see text. 



Floquet side-bands in the LDOS structure with respect 
to an enhancement of the hopping amplitude T. Dis- 
tinct gap-states evolve which induce a transition from the 
Mott-insulator state to the liquid phase. Both features 
result in severe changes of the fermionic band structure 
and therefore cause significant changes of e.g. optical and 
conduction properties. 

The behavior of the density of states as function of the 
external modulation frequency exhibits two limiting 
cases with a cross-over regime between them. For the 
limit of small modulation frequencies VLl ^ it is in- 
teresting to note that all Floquet-modes gain more and 
more spectral weight. This signals the onset of an or- 
thogonality catastrophe, as predicted by P. W. Anderson 
(AOC) [325. The system undergoes a transition to a new 
ground state which is orthogonal to the original equilib- 
rium state, caused by the change in the potential of the 
system Fig. O Technically this is seen as drastic enhance- 
ment of the arguments of the Bessel functions, e.g. for 
the non-interacting expression in Eq. (jlj. Considering 
the behavior of the Bessel functions, we note that with 
decreasing modulation frequency the number of Floquet 
modes, which contribute, increases. Numerical calcula- 
tions however are always limited to treatments regarding 
a finite number of Floquet modes. Therefore our numer- 
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ics is optimized towards the limit of maximum validity for 
an affordable amount of computational effort. A detailed 
analysis can be found in Fig. [3] and caption. 

In Fig. [5] we discuss the behavior of the LDOS and 
the occupation number for increasing external modula- 
tion frequency Ql- For frequencies Ql < ^ the behavior 
of the ultracold Fermi gas changes from Mott-insulating 
to Hquid or conducting. Pronounced Floquet-side bands 
[26| are developed and intersect in between the Hubbard- 
bands. The Mott-gap almost disappears and a liquid or 
conducting regime is established, where the liquid den- 
sity of states can be continuously driven by the external 
modulation. Regarding the occupation number of the 
gap states we find a step-like behavior for long wave- 
length modulations which can be interpreted as the ab- 
sorption or emission of energy quanta (phonons). In the 
non-equilibrium fermionic distribution function we de- 
rive that for long wavelength modulations the majority of 
fermions resides in states below the Fermi edge (a; = 0). 

The two limiting regimes of small and large lattice 
modulations are separated by a cross-over at — 1 
(see Fig. [2]). At the crossing, the modulation induced 
Floquet-side bands are forced to intersect in the gap and 
acquire a maximum of spectral weight (see Fig. H]) in 
this area. Moreover the occupancy from the states right 
above the lower Hubbard-band is shifted towards states 
right below the upper Hubbard-band and additionally 
the entire gap is almost equally occupied. We further 
remark that the excitation behavior to reach the upper 
Hubbard-band at the crossing changes from virtual, i.e. 
successive absorption, to direct. 

Right above the crossing, for external modulation fre- 
quencies > 1 we find, that the liquid behavior is dra- 
matically reduced. In the lower panels of Fig. [5] we dis- 
cuss the LDOS and the distribution function in the for 
external frequencies above = 1- We have to remark 
that above the crossing the LDOS shows a significant 
change in the gap. The step-like structure vanishes and 
the system reveals for further increasing almost the 
behavior of a Mott-insulator with an unconventional oc- 
cupation number in the Mott-gap associated with a weak 
spectral weight there. Combining both, the discussion 
concerning the spectral weight and the distribution num- 
ber right above = 1 we find a pronounced shift of oc- 
cupation which resembles a "trapping" of particles right 
below the upper Hubbard-band, which is an "inversion" 
alike situation. The Fermions can not relaxate because 
no unoccupied states are within reach if we assume only 
emission processes of an integer number of phonons. For 
significantly faster lattice vibrations of the confining po- 
tential the system is not able to follow the perturbations 
and returns to a stationary state similar to equilibrium. 
The trapping of the occupation is preserved. 

To conclude we discussed a theory of ultracold 
fermionic atoms described by a Hubbard model includ- 
ing strong repulsive interactions. Periodic lattice poten- 



tial modulations which drive this system out of thermo- 
dynamical equilibrium were considered. By investigat- 
ing a Floquet-Keldysh Green's approach we find a cross- 
over between two limiting characteristics, the AOC for 
l^L ^ and a quasi-equilibrium solution for oo. 
Pronounced side-bands lead to a rather complicated den- 
sity of states in the gap for < 1 which indicates a tran- 
sition to the Hquid or conducting regime right at the onset 
of the modulation. In the vicinity of the cross-over we 
find a maximum of spectral weight inside the original ex- 
citation gap. For external frequencies Ql > ^ population 
trapping in the gap is observed. Beyond the system ap- 
proaches an equilibrium-like Mott-insulator regime which 
indicates that the system is not able to follow very fast 
perturbations. 
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